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Abstract 

The closed analytical expression for the Uehling potential is derived. The Uehling potential 
describes the lowest-order correction on vacuum polarization in atomic and muon-atomic systems. 
By using our three-term analytical formula for the Uehling potential we evaluated the vacuum 
polarization correction(s) for the ground states in some two-electron ions. The analytical formula 
for the Uehling potential allows us to derive the general expression for the interaction potential be- 
tween two point electric charges. This analytical formula provides the correct asymtotic behaviour 
at large distances between two interacting electric charges. 
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I. INTRODUCTION 



Various effects related to the vacuum polarization in atomic systems have been analyzed 
by many famous theorists, including Dirac, Heisenberg, Pauli and others. A large number 
of fundamental papers of these authors were stumulated by investigations of the vacuum 
polarization problem (see, e.g. ^ and j^). Briefly, the vacuum polarization in atomic 
systems can be described as the creation of virtual electron-positron pairs by a photon 
moving in the field created by an electric charge eQ (see, e.g., {4], js]). These virtual 
(e~, e"^)— pairs will interact with the original photon, and such an interaction will affect 
its propagation, an effect termed the vacuum polarization by Uehling in 1935 {g]. It can 
be detected in all atoms and muonic atoms. For heavy atomic systems, e.g., for the Pb, 
Bi and U atoms and corresponding muonic atoms, the effect of vacuum polarization has 
been discussed in various papers (see, e.g., 

a, a 

and references therein). There are some 
advantages in studying the effect of vacuum polarization in few-electron atoms and muonic 
atoms. Indeed, in the case of light atomic systems with relatively small nuclear charges 
Q < 10, we can assume that in the vicinity of an atomic nucleus with positive charge eQ 
the actual interaction potential will slightly be different from the scalar Coulomb potential 
created by the two point particles. It is found in this study that such a potential can 
still be written in a closed analytical form. In reality, this means that the pure Coulomb 
potential must be replaced by some 'more realistic' potential such as V{r) = ^^{1 + f{Q, r)), 
where f{Q, r) is the function which represents the deviation of the actual potential from the 
Coulomb potential at small r. Note that this approach (Dirac's approach) corresponds to the 
non-relativistic picture, but it can easily be generalized to the relativistic case. In general, 
the function f{Q,r) is determined with the use of methods of Quantum Electrodynamics. 
The main goal of this study is to derive the closed analytical expression for the f{Q, r) 
function. 

The expression for the f{Q,r) function can be derived explicitly, if we know the low- 
est order radiative corrections to the photon Green's function and for the vertex function. 
Formulas for radiative corrections to the photon Green's function are very well known (see, 
e.g., jsl). By choosing the lowest order correction one obtains the modified Coulomb law 
in vacuum. This problem was considered for the first time by Uehling in 1935 [6] and the 
associated f{Q, r) function which arises in this case is called the Uehling potential. Since 
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early studies of vacuum polarization many different people tried to derive the closed ana- 
lytical expression for the Uehling potential. Unfortunately, such an expression has not been 
found and it substantially delayed theoretical research of the vacuum polarization. In this 
study we have derived the closed, three-term analytical formula for the Uehling potential. 
This formula is used to obtain the properties of the Uehling potential and perform numerical 
evaluations of the corresponding corrections to the total energies of bound states in light 
two-electron atoms and ions. The closed analytical formula for the Uehling potential is the 
main result of this work. Note that our formula includes only three terms each of which con- 
tains either the modified Bessel function of the zero order KQ{2br), or the recursive integrals 
of this function, i.e. the Kii(2br) and Ki2{2br) functions. This formula provides numerous 
advantages for theoretical analysis as well as for numerical computations. 

Note, however, that the Uehling potential has wrong long-range asymptotics, i.e. its 
behaviour at large r is not correct. Therefore, it must be modify to represent the interac- 
tion between the two point electric charges correctly at arbitrary, in principle, interparticle 
distances. This problem has substantial scientific and methodological interest. Indeed, in 
classical electrodynamics the interaction V between two point electric charges qie and q2e 

2 

is described by the Coulomb law, i.e. V ~ ^^f^- This expression for V does not change 
its form in the non-relativistic quantum mechanics [9i] , when the Planck constant h is finite. 
However, if we also assume that the speed of light c is finite, then Quantum Electrodynam- 
ics leads us to the necessaty to modify the Coulomb law. The main correction is related 
to the vacuum polarization in the spatial areas close to the electric charges. In the lowest 
order approximation such a correction is represented by the Uehling potential js], which is 
correct at short interparticle distances. However, the Uehling potential does not provide 
the correct long-range asymptotics, i.e. its behaviour at large r is not correct. Therefore, it 
must be modify to represent the interaction between the two point electric charges correctly 
at arbitrary, in principle, interparticle distances. This problem has substantial scientific 
and methodological interest. In particular, the Coulomb potential modified by a number 
additional terms, which represent small QED-corrections, can be used directly in the highly 
accurate calculations of few-electron atoms and ions. On the other hand, it is interesting to 
see the transformation of the usual Coulomb potential under the rules of modern Quantum 
Electrodynamics. 

The explicit derivation of such a corrected Coulomb potential modified by a number 
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additional terms is the second goal of this study. Such a potential can directly be used 
in highly accurate calculations of few-electron atoms and ions. On the other hand, it is 
interesting to see the transformation of the usual Coulomb potential under the rules of 
modern Quantum Electrodynamics. 



II. UEHLING POTENTIAL 



As mentioned above, in the lowest order approximation the vacuum polarization is de- 
scribed by the Uehling potential jo] which is represented as a small, short-range correction 
to the regular Coulomb potential acting between the two point particles which have electric 
charges e and Qe. Formally, the contribution from such a correction can be represented 
as aQ{r), where a = is the fine structure constant and Q{r) is the regular Coulomb 
potential acting between two charged particles. The explicit formula for the Uehling po- 
tential is derived by computing the lowest-order radiative corrections to the photon Green's 
function. In fact, if the current Ju{x) is known (here and below ly = 0, 1, 2, 3), then the 
electromagnetic field Fy{x) which is created by this current is obtained from the following 
formula 

F,y{x) = z J G^^J{x — x')J^{x')dx' (1) 

where G^^J(x) is the photon Green's function which also contains the lowest order radiative 
correction. The integration in Eq.(IT]) is performed over all 4-dimensional Minkovskii space 
where we have g'^^ = (1, —1, —1, —1) and Qfj^u = 0, if ^ u. 

If the field is created by a point electric charge which is at rest, then Ji, = teQSiyoSir) and 
the corresponding expression for scalar potential 0(r) takes the form 

0(r) = leQ r G^^\x)dt = / G^^\\i, 0) exp(2k ■ Y)(f\<i (2) 

J-oa J 

where x = {ct, r) and the last integral is a three-dimensional integral over k-space. Here we 
designate 0(r) = Jo{r) = Joir), where r =| r |. The analytical formula for the G*-'^''(k, 0) 
function is well known (see, e.g., jsj and references therein) 

G(^)(k,0) = + — [\(l-x)-\n\l + -^x(l-x)]dx] (3) 
^ n Jo L m"' J J 

2 

where a = is the fine structure constant. In the relativistic units with c = 1 and 
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h = l,a = By substituting this expression in Eq.([2]) one finds the formula 



eQ Qea v^{l 



Anr 7r(27r)'^ jo 



V 



exp(ik ■ r) 

1— ti'^ 



where = 2x — 1. The integral over the k variable in this formula is 



exp(zk ■ r) 

^2 _|_ 4m2 

1 — 1)2 



27r2 



exp(- 



2mr 



(4) 



(5) 



Therefore, the potential (p{r) in Eq.(j4]) can be written in the form 



exp 



2mr 



dv 



(6) 



where we have used the fact that the expression under the integral is an even function of v. 



The integral in Eq.(l6l) can be transformed by introducing the new variable ^ 
From the definition of this variable ^ one finds 



1 



vT—v- 



and 
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.(1 + ^ 
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By using these formulas we reduce the integral from Eq.([6]) 



to the form 
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where a = 1^ is the fine structure constant (in relativistic units). The sum of Eq. (fTOl) with 
the Coulomb potential ^ is 
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Anr 
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(11) 



The energy of an elementary charge e in the field of this potential is written in the form 



e(j){r) = a 



Q 



1 + 



2a 



3tt Ji 



exp(-2m^r)(l + 



1 xv/F^ 



dC 



(12) 



In the case of atomic systems, where m = rrie, one finds the following expression for the 
ec(p{r) interaction energy (in atomic units h = 1, rrif. = 1 and e = 1) 

Q 



ec0[r] 



Q 
r 



1 + 



2a 

37^ 



+ 00 



exp{—2a ^^r)[ \ + 
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1 ^ v^e^ 
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+ Uir) 



(13) 



where Q is the integer product of two particle charges, c0 is the 0-component of the 4- 
dimensional vector-potnetial = (c0, A) and a is the fine structure constant {a^^ ^ 
137.03599911, see, e.g., [10]). The 'additional' potential energy U{r) in Eq. (fT^ corresponds 
to the contribution from the so-called Uehling potential which is generated by the lowest or- 
der radiative correction to the electrodynamic Green's function, or photon Green's function, 
for short. As follows from Eq.( !T3|) the Uehling potential differs from zero only for very short 
distances r < Ae, where Ae = = aa^ ^ 3.861592 ■ 10~^^ cm is the Compton wavelength 
of the electron and oq is the Bohr radius. 

In our derivation of analytical formulas for the U{r) term in Eq. flT^ all important de- 
tails are shown explicitly. The formulas given above allow one to understand some mis- 
takes/misprints which are often made in the consideration of the Uehling potential. For 
instance, the integral representation (Laplace tranform) of the Uehling potential given in 
0] contains an extra negative power of the ^—variable. However, our main goal below is to 
derive and discuss the closed analytical expression for the Uehling potential. According to 
the modern scientific literature such a formula does not exist and only the corresponding 
integral transformations are used in actual applications. 



A. Long-range and short-range asymptotics of the Uehling potential 

As is well known (see, e.g., js] and references therein) the integral 

r"^ 1 \ \/P - 1 

I = J^ expi-2m^r){l + —)^^^^d^ (14) 

can be computed analytically for very small and very large values of its parameter mr. For 

mr <^ 1 the approximate value of this integral is 

5 

I ~ — - — 7 — In(mr) (15) 

or, for m = me and in atomic units 

/~---7 + lna-lnr^ -6.33079265675 - In r (16) 
6 



where 7 ^ 0.5772156649 ... is the Euler constant (see, e.g., [ll|) and Ina ^ -4.92024365857. 

In the opposite case, i.e. when mr ^ 1 the long-range asymptotics of the integral /, 
Eq.dUD, is 

3v^ exp(-2mr) 
^ - ^ : — 7F— (17) 



or, for m = rup and in atomic units 



S-v/tt 3 exp(— 2a! ^r) . exp(— 2a ^r) , . 

■ a5 ■ 3 ^ 4.1433745897 ■ 10"^ ■ (18) 

8 /■ 2 f2 



where a^^ ^ 137.03599911 (see, e.g., 



10|). 



By using these formulas for /(r) we can find the appropriate analytical expressions for 
the long-range and short-range asymptotics of the Uehling potential U (r) and its sum with 
the Coulomb potential. Indeed, the Uehling potential is U{r) = ^^I{r), while for the total 
interaction potential we have 0(r) = ^ 1 + ^I(r) . Here and everywhere below in this 
work we shall apply only atomic units. From these expressions one can find the explicit 
short-range and long-range asymptotics of the U{r) and 0(r) potentials. Note that these 
asymptotic formulas for the U{r) potential have been known since the middle of 1930's 6| 
(see discussion in [l8|). A number of people tried to produce the general formula for the 
Uehling potential. Such a formula must provide the correct asymptotics at latrge r and small 
r. All such attempts were unsuccesfuU and around the middle of 1950's it was finally assumed 
that simple analytical expression for the Uehling potential U{r) at arbitrary interparticle 
distances r does not exist. Since then only numerical methods are used for actual atomic 
and muon atomic systems. In the next Section we derive the explicit analytical formula for 
the Uehling potential U (r) for all values of r. 



III. ANALYTICAL FORMULA FOR THE UEHLING POTENTIAL 

In numerous books and textbooks on Quantum Electrodynamics (see, e.g., jsl, [l3|) one 
can find a statement that the integral in Eq. (|T3l) 'cannot be evaluated in closed form but is 
easily solved numerically'. In other words, it is widely assumed that the analytical expression 
for the U{r) potential does not exist. In order to show that this statement is incorrect we 
proceed to obtain the closed analytical form for the Uehling potential U{r). First, consider 
the following integral 

2^2 J ^2 

where a is some positive numerical parameter. By using the substitution ^ = cosh a; we 
reduce lu to the form 

1 1 \ . . . 1 . . 1 



r+'^ / 1 \ - 1 

^u{a) = l exp{-aO{l + ^)^^^^d^ , (19) 



ea;p(-acoshx)(^l - ^ ^^^^^ ^ ~ 2cosh'^x ^'^^ ^ ^^^^^ " 2'^^^*^"'* ~ 2^^ii'^) 
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where KQ{a) is the modified Bessel function of zero order, i.e. 



^2k 



POO 

Ko{z)= exp{-zcosht)dt = Y.{tlj{k + l) + ln2-\nz) , 
Jo ^0 22'=(A;!)2 

where is the Euler psz-function defined by Eq.8.362 from \v^. The functions Ki2{z) 
and Kiti{z) in Eq. (l2U]) are the recursive integrals of the Kq{z) function, i.e. 

poo roo 

Kii{z) = I KiQ{z)dz , and Kin{z) = / Kin-i{z)dz , (21) 

J z J z 

where n >1 and Kio{z) = Ko{z). By using the relation given by Eq.11.2.14 from IqI for 
the recursive integrals Kin{z) (where n > 0) one can reduce Eq. fl20l) to the following form 

Iu{z) = {l+'^)K,{z)-^K^,{z) + {\-'^)K^,{z) . (22) 

In the present case we have z = a = 2br, and therefore one finds the following explicit 
formula for Uehling potential U{r) = U (2br) in position space 

U{2br) = ^[(l + ^^)M2br) - jKt^{2br) + (-^ + l)Kt,{2br)] , (23) 

where the notation Q stands for the electric charge of the nucleus and b = a~^. Now, Eq. (fT3|) 
takes the form 



^ + U{r) (24) 
r 



This is the general form of the Coulomb interaction energy between electric charges Qe and e 
which includes the lowest-order vacuum polarization correction. As expected this expression 
depends only upon the fine structure constant, i.e. upon the h and c universal constants and 
electric charge e (and Qe), but it does nor include any particle's mass. The generalization of 
these formulas to the case of two interacting electric charges q\ and is obvious and simple. 
Note also that the two limiting forms of the Uehling potential mentioned above (see. Section 
II) which correspond to the cases when r ^ aa^ = Ag and r 3> aa^ = Ag, respectively, can 
be obtained from these expressions Eqs.(l23l) - IHM . 

Note that the formula Eq. (l23l) includes only three different terms each of which contain 
either the modified Bessel function of zero order KQ{2br), or one of the two recursive inte- 
grals from this function, i.e. the Kii[2br) and/or Ki2{2br) functions. The formula, Eq.( l23|) . 
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has many advantages for theoretical analysis as well as for numerical computations. Indeed, 
the parameter b in this formula equals q;^^(^ 1), and therefore, Eq. fj23|) can be used to 
obtain various approximate formulas for the Uehling potential. Each term in such formulas 
has a different order upon the fine structure constant a (or a^^). Moreover, this closed 
(and relatively simple) analytical formula is useful in many other problems, e.g., to perform 
secondary quantization of Maxwell equations which already include vacuum polarization cor- 
rection. In numerical computations the formula, Eq.f l23p . allows one to avoid very painful 
and time-consuming procedures of numerical integration. Briefly, this means that each ma- 
trix elements of the Uehling potential is computed analytically in an arbitrary, in principle, 
basis of radial few-body functions. The closed analytical formulas cannot lead to the loss of 
numerical accuracy in computations. Furthermore, in this case we do not need to investigate 
the convergence of the matrix elements and overall results upon various parameters arising 
in the methods used for numerical integration. 



IV. RESULTS AND DISCUSSION 

To check our analytical formula, Eq.f l23p . for the Uehling potential we determine the cor- 
rections to the bound state energies of some atoms and ions associated with such a potential. 
In other words, we determine the lowest order vacuum polarization correction(s) for some 
simple (two-electron) atoms and ions. In general, to evaluate the correction to the bound 
state energy associated with the Uehling potential one needs to solve the non-relativistic 
Schrodinger equation for an atomic system with the pure Coulomb potential to very high 
accuracy. At the second stage of the procedure the non-relativistic wave functions are used 
to compute the expectation values of the Uehling potential. Formally, in this approach we 
have to consider the corrections which correspond to all levels of the perturbation theory, 
e.g., the lowest-order corrections, second-order corrections, etc. The convergence of such 
series from perturbation theory is often a very complex question. An alternative method 
is based on the direct solution of the non-relativistic Schrodinger equation with the new 
interaction potential which includes the Uehling potential. The corresponding correction 
is determined as a difference of the two energies (eigenvalues) obtained from the solution 
of the two Schrodinger equations: one of these equations contains the pure Coulomb inter- 
action potential, while in the second equation the interaction potential is represented as a 
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sum of Coulomb and Uehling potentials. The wave function obtained in the result of the 
numerical solution of the Schrodinger equation already contains all corrections related with 
the vacuum polarization correction. By using this method below we determine the vacuum 
polarization corrections (or Uehling correction) for the two-electron helium atoms and for 
some helium-like ions. 

In our calculations we have used the highly accurate variational wave functions con- 
structed from the exponential basis funcions of three relative coordinates r32,r3i and r2i. 
Each of these trial wave function includes 400 terms (exponents) with the carefully opti- 
mized non-linear parameters (more detail about exponential basis functions can be found, 



e.g., m 



15| and references therein). Finally, these trial wave functions for each two-electron 



atomic system have very good numerical accuracy. It provides ^12-14 stable decimal 
digits in the total energy. It is clear that such an accuracy is sufficient to evaluate the lowest 
order vacuum polarization correction directly and to high accuracy. The results of numerical 
evaluations of the vacuum polarization correction for the °°He, ^He and ^He atoms can be 
found in Table I. Our results in Table I include the non-relativistic Coulomb energies Enr 
computed with = 400 basis wave functions and = 4000 basis functions (exponents) 
[igI . The last values are designated in Table I as i?„r(oo), since it is very close to the exact 
answer. Table I also includes the total vacuum polarization correction AEueh for all three 
interparticle (Coulomb) interactions in the He atoms. We also determine the vacuum po- 
larization corrections for the electron-nucleus and electron-electron interactions separately. 
They are designated in Table I as Ai?^^ and /S.E^^^, respectively. 

As follows from the results of Table I, the overall correction to the energy due to inclusion 
of the vacuum polarization correction for the two-electron He atom is relatively small ~ - 
7.2555-10"'' a.u. However, such a difference can easily be detected in modern high-precision 
optical experiments. The contribution of the vacuum polarization correction for the electron- 
electron interaction is substantially smaller -1.1051-10"^ a.u. In other words, in the first 
approximation we can neglect the correction for the electron-electron interactions in all few- 
electron atoms. Another interesting fact, which follows from Table I, is the mass dependence 
of the corrections. The infiuence of variations in the masses of the helium atom(s) °°He, "^He 
and ^He appears in the 4th decimal digit in the vacuum polarization correction. 

The results of numerical computations of the lowest order vacuum polarization correc- 
tions in some two-electron ions can be found in Table II. In Table II we consider the ground 
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1^S{L = 0)— states of the °°Li+, ""Be^"*", "^B^"*" and "^C^"*" ions. Vacuum polarization correc- 
tions AEueh have been determined for each of these ions and shown in Table II in atomic 
units. As follows from Table II the vacuum polarization correction rapidly increases with 
the nuclear charge Q. Note that our results for the lowest order vacuum polarization correc- 



br the ground state of the He atom concide well with the results of earlier 



tion(s) computed 

studies (see, e.g., [17| and references therein). In the future work we are planning to apply 
the formula, Eq.f l23p . to few-electron atoms and ions. The matrix elements of the Uehling 
potential must be obtained in different systems of radial basis functions known for few-body 
systems. Another interesting problem is to compare the efficiency of our method based on 
the analytical formula, Eq.f l23p . for the Uehling potential with the methods based on the 
numerical integration. In particular, such a comparison can be made for the two-electron 
atomic systems discussed above. 

V. CORRECT ASYMPTOTICS AT LARGE DISTANCES 

Our analytical formula for the interaction potential between the electron and atomic 
nucleus (with the electric charge Qe) is written in the following form (in atomic units h = 
1, nie = 1 and e = 1) 

= 9. + u{r) (25) 
r 

where U{r) is the Uehling potential, Kq{z) is the modified Bessel function of zero order, 
while Kii{z) and Ki2{z) are the recursive integrals of the Kq{z) function (see above). Note 
that the formula, Eq. (|25l) . is corret mathematically, but it has a wrong 'physical' asymptotic 
at large distances r, e.g., at r > lOaao, where = ^^^-^ is the Bohr radius. First, note that 
the short range asymptotic of the (/)(r) potential, Eq.f l24p . takes the form 

0(r) 1,^0=^ + [-^ - 27 + 2 In a - 2 In rl I (26) 



where 7 ^ 0.5772156649 ... is the Euler constant (see, e.g., (llj) and Ina ^ -4.92024365857, 



The long-range asymptotics of the potential 0(r), Eq. (l2^ . is (in atomic units) 

«{l + ^exp(-2.,} (27) 
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This means that the long-range asymptotics of 0(r) decreases with r exponentially. This is 
not correct physically, since the exponential function vanishes at r — t- +oo very rapidly. In 
reality, such an incorrect asymptotics of the interaction potential 0(r) at large distances r 
is overweighted by the lowest order QED correction to the electromagnetic field (E,H), or 
to the pure electric field E in our case, when H = 0. This correction provides the correct 
power-type dependence of the interparticle potential at large disctances. It is directly related 
to the non-linearity of the Maxwell equation for the EM-field 

By following the method used in 18|] let us obtain the lowest order QED correction to the 
pure electric field E. Consider one point electric charge located at the center of coordinates 
X = 0,y = 0, z = 0. This charge is assumed to be at rest. The corresponding Lagrangian L 
which includes the lowest order QED correction is written in the form (in regular units) 

where E = — V(;/)(r). By varying the potential 0(r) in the equation (the fundamental equation 
of least action) (jj] 

6 J L[(f){r)ydr = (29) 

one finds the following differential equation 

d(l) ^ e^h .dc/yy _ C 
dr QOTT^mgC'' ^dr ' 

where C = —j^- Assuming that (/)(r) = ^ + ^jj{r), where ilj{r) is a very small correction to 
the Coulomb potnetial, we reduce the last equation to the form 



dr gOvr^mfc^ ^64n^r^ 

From this equation one finds 



Therefore, the total interaction potential 0(r) is 



V^W = -T,7;5:i3(^) (32) 



Qe e^Ti 
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225™fc7^647r3/r4 
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In the relativistic units ^ = 1, c = 1 and a = ^ the last formula takes the form 

' in 



(33) 
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(34) 
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This expression is the sum of the Coulomb potential and small correction which exactly 



coincides with the well known Wichmann-KroU potential 
expression is 



18| . In atomic units the last 



r = 

r 



1 - 



2257rr4 

and Wichmann-KroU potential Wk{t) is 



3^,7 



r 2257Tr^ 



(35) 



Hr) = W,ir) = (36) 

Note that this potential is not regular at the origin. However, as follows from the Appendix 
such a singularity is formal, since Eq. fl5U]l is not correct at r — )■ (see Appendix). It 
can be removed by using the substitution r — )■ r + a in atomic units or r — r + aao in 
regular units. This substitution removes singularity of the Wxir) potential at r = and 
changes its behaviour at short interparticle distances r < a. However, at such distances 
the Wichmann-KroU correction is not important and its contribution is significantly smaller 
than contribution from the Uehling potential. On the other hand, at large distances, i.e. 
for r > lOaao (and even for r > aao) the contribution from the Wichmann-KroU potential 
plays the leading role. Finally, we can write the correct expression for the Wichmann-KroU 
potential Wxir) 

Note that this potential is always negative. In general, at r ^ oq this potential is 
very small by its absolute value. However, its overall contribution rapidly increases with 
the nuclear charge Q. Furthermore, it decreases with the distance r as ~ r~^, i.e. non- 
exponentially. It is clear that at some large distance the Wichmann-KroU correction will 
always exceed the contribution form the Uehling potential U (r), Eq. fl2^ . Thus, the following 
interaction potential (in atomic units) of the two point electric charges 

— + U{r) + WKir) = — + 

r r Snr 



$(,) = « + U(r) + W,ir) = ^ + . [(1 + :fl)A-o(^) - f /«.(^^ 



+ (- - ^)K^2(-)] - "'f ^. (38) 
Vfi .Sn-sy ^\ryJ\ 2257r(r + a)5 ^ ' 



'5 r2 X .2rNi 2QW 
3a;^^ ^ a 

has the correct asymptotic behaviour both at small and large interparticle distances. In 
the case of interaction between two point electric charges gic and we need to replace in 
Eq.( l38i) the factor Q by the product qiq2- The potential $(ri2) takes the form 



$ ^2 = + ^ [l + —]Ko[ ]-—Kii[ 

ri2 onr 12 3q;^^ ^ a ^ ba ^ 
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V6 Sa^y a 2257r(ri2 + a)^ ^ ^ 

where ri2 is the distance between particles 1 and 2. 

Thus, we have derived the formula for the interaction between two point electric charges 
which includes the lowest order QED corrections to the Coulomb potential. Our formulas, 
Eqs.f l38p - (!39|) . provide the correct asymptotic behaviour at arbitrary interparticle distances. 
The formulas, Eqs. (l38l) - ( l39l) . can directly be used in highly accurate computations of the 
bound states in few-electron atoms and ions. The explicit computation of the electric field 
strength E, i.e. the spatial gradient from the $(ri2) potential, Eq.( l39|) . is straightforward, 
but the Fourier resolution of the corresponding 'electrostatic' field is extremely complicated. 
Note also that for the potential $(r), Eq. (!38|) . the condition A$(r) = 47rQ^(r), where 
A = V ■ V = div{grad...), is not obeyed. 

VI. CONCLUSION 

Thus, we have derived the closed analytical expression for the Uehling potential which 
represents the lowest order vacuum polarization correction(s) in atomic systems. It is shown 
that the Uehling potential can be represented as a sum of the modified Bessel function of 
the zero order Ko{2br) and the recursive integrals of this function, i.e. the Kii{2br) and 
Ki2{2br) functions (see Eq. fl23|l ). Based on our explicit and relatively simple three-term 
formula, Eq.f l23p . we can reject a well known existing conclusion that the Uehling potential 
cannot be represented by the closed analytical formula. Moreover, we have found that such 
a simple analytical formula drastically simlifies theoretical analysis of vacuum polarization 
in light atoms and ions. Numerical algorithms based on this formula are usually very fast 
and effective for accurate calculations of atomic systems. Briefly, we can say that the use of 
this closed analytical formula for the Uehling potential opens a new avenue in investigation 
of the vacuum polarization in light and heavy atoms and ions and muonic atoms and ions. 
It can be used for various theoretical purposes, e.g., to perform secondary quantization of 
Maxwell equations modified to the form which includes the lowest order vacuum polarization 
correction. 

To check our analytical formula for the Uehling potential, Eq. (!23|) . we have determined 
the lowest order vacuum polarization correction for some atoms and ions. In particular. 
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these corrections are evaluated numerically for the ground l^S"— states of the two-electron 
He atoms and He-like ions. The overall contribution from vacuum polarization for He-atoms 
is evaluated as ~ -7.2555 -10^^ a.u. The main contribution to this correction comes from 
electron-nucleus interactions. The effect of vacuum polarization for the electron-electron 
interactions is very small (at least by an order of magnitude smaller) and can be neglected 
in the lowest-order approximation. We also determine the vacuum polarization correction for 
some two-electron (or helium-like) ions. In particular, we have determined such corrections 
for the ground l^iS— states in the Li+, Be^'^, B'^^ and C^^ ions. 

We also derived the general formula, Eq.f lS^ . for the interaction between two point 
electric charges. This expression includes the lowest order QED corrections to the regular 
Coulomb potential. Our formulas, Eqs.f lSK]) - fl5^ . have the correct asymptotic behaviour at 
small and large interparticle distances. Briefly, we can consider this formula as the modified 
Coulomb potential which agrees with the rules of modern QED in the lowest order. 
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Appendix 

The equation, Eq. (!30|) . is a qubic equation in respect to the radial derivative of the 
unknown potential '0(r). It can be re-written in the form 

^y' + y + ^.=^ (40) 



where U = As follows from the general theory of qubic functions this equation has only 
one real root (the discriminant of this equation is negative). In fact, this equation is easily 



reduced to the form of a monic trinomial (see, e.g., [19|) 

l/'+Pl/+^ = (41) 

where p = — and q = 



The Cardano method gives the only real root of Eq. (jH 



q I q'^ 3 



q I q^ 



\ 2r2 V 4r4 27 ^ 2r^ V 4r4 27 
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+ ^ (42) 



In other words, we obtain the following differential equation 

d(f) 3 
dr \ 

At r — 7- oo this equation has the solution 'tp{r) = 2(|)5 ■ r + B, where B is some numerical 
constant. Sucn an asymptotics is not correct. At r — the analogous solution takes the 
form ijj(r) = —3q^ ■ + C, where C is some numerical constant. The numerical constant 
C must contain the fine structure constant a. The overall contribution of the Wichmann- 
Kroll term at small r must be smaller than the analogous contribution from the Uehling 
potential. This gives C ^ a", where n > 2. Now, it is easy to check that the substitution 
r—j-r + a (or r—)'r + aao) used in the main text is indeed the correct way to regularize 
the Wichmann-KroU potential at small r. 
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TABLE I: The total non-relativistic energies £;„r(400) of the ground 1^S{L = 0)-state of the °°He, 
^He and ^He atom and vacuum polarization corrections (in atomic units). 



atom 


°^He 


^He 


^He 


Enr (oo) 


-2.9037243770341195831 


-2.9033045577295687858 


-2.9031672107105772470 


^nr (400) 


-2.90372437703405 


-2.90330455772950 


-2.90316721071051 




-7.255456205-10-^ 


-7.252434438-10-^ 


-7.251448021-10-^ 




-7.366031256-10-^ 


-7.362968355-10-^ 


-7.361966472-10-^ 




1.10575051-10-*^ 


1.10533917-10-^ 


1.10518451-10-^ 



The energies AEf^'^f^ include the vacuum polarization correction for the electron-nucleus 
interactions only, i.e. AE^^ = AEueh - AE^f^. 

TABLE II: The total non-relativistic energies £;„^(400) of the ground 1^S{L = 0)-state of the 
°°Li+, °°Be^+, °°B^+, °°C^+, and vacuum polarization corrections (in atomic units). 



ion °°Li+ ~Be2+ °°B^+ ~C''+ 

Enr{oo) -7.27991341266930596 -13.65556623842358670 -22.03097158024278154 -36.40624660189853031 
E„^(400) -7.27991341266315 -13.6555662384085 -22.0309715801925 -36.4062466018083 
AEueh -2.706070460-10-6 -6.714805549-10-^ -1.342648935-10-^ -2.349488675-10-^ 
AE^'^f^ -2.762425832-10-6 -6.873080812-10-^ -1.377060071-10-^ -2.413281437-10-^ 
AE^lf^ 5.6355372-10-^ 1.5927526-10-^ 3.4412136-10-^ 6.3792762-10-^ 

The energies AE^^^ include the vacuum polarization correction for the electron-nucleus 
interactions only, i.e. AE^^ = AE^eh - AE^lf^. 
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